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RUNAWAY  ELECTRONS  IN  COLLECTIVE  ELECTRIC  FIELDS  ; 

i 


I INTRODUCTION 


Recently,  there  have  been  numerous  reports’  of  runaway  electrons  in  tokamaks.  The  energy 
of  runaway  electrons  observed  in  low  density  discharges  (n  < 10 cm  is  lower  than  that  observed 
in  higher  density  discharges,  it  has  also  been  observed  that  the  synchrotron  radiation  power  often  in- 
creases in  steps,  and  is  accompanied  by  radiation  near  and  below  the  plasma  frequency  As  first 

pointed  out  by  Coppi  et.  al.,’  the  reduced  plasma  frequency  radiation  (tu*  =k,|//c<Opp)  in  Alcator 
could  be  explained  by  an  anomalous  Doppler  shifted  mode,^  where  is  the  wave  frequency,  k is  the 
wave  number,  and  k||  is  its  component  in  the  toroidal  magnetic  field  direction.  Pitch-angle  scattering 
of  the  runaway  electrons  by  the  anomalous-Doppler-shifted  mode  can  lead  to  an  increase  of  synchro- 
tron radiation. Papadopoulos  et  al.'^  also  show  that  pitch-angle  scattering  can  cause  a bump  to  form 
on  the  tail  of  the  electron  distribution.  The  bump-on-tail  instability  will  then  be  excited  and  give  rise 
to  the  Wp,.  mode  which  is  observed.  Moreover,  anomalous-Doppler-shifted  and  bump-on-tail  modes 
may  cause  many  runaway  electrons  to  become  trapped  in  local  mirrors,"  VS  drift  to  the  wall  of  the 
vacuum  vessel,  and  destroy  its  liner.  Such  liner  destruction  may  have  been  observed  by  the  TFR 
group. ^ The  macroscopic  behavior  of  the  anamolous-Doppler-shifted  mode  is  described  in  detail  by 
Choi  and  Horton. " 


I 

I 


In  this  paper,  we  study  the  interactions  between  the  runaway  electrons  and  the  collective  modes 

using  a Fokker-Planck  model  including  a self-consistent  quasilinear  diffusion  term.  First,  we  show  that 

if  only  the  anomalous-Doppler-shifted  mode  is  present,  it  gives  rise  to  the  bump-on-tail  distribution. 

This  is  because  of  a velocity  threshold  for  the  anomalous-Doppler-shifted  mode.  Second,  we  show  that 

the  combination  of  the  anomalous-Doppler-shifted  mode,  the  bump-on-tail  mode  and  a weak  toroidal 

electric  field  Fn  , can  lead  to  synchrotron  radiation  similar  to  that  experimentally  observed.  In  section 

If,  we  review  the  analytical  theory.  In  section  III,  we  present  the  numerical  model.  The  results  arc 

given  in  section  IV,  and  these  results  are  summarized  in  section  V. 

Note:  Manuicript  lubmitted  September  14,  1977. 
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II  ANALYTICAL  MODEL 


For  electrostatic  waves  in  a homogeneous  magnetic  field  with  immobile  ions,  we  use  the  Harris 
dispersion  relation: 


* eo^  I dv  y j 

•>  1 Acii  V|i  - 
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(I) 


where  is  the  Bessel  function  of  the  first  kind  of  order!,  and  w,.,.  is  the  electron  cyclotron  frequency. 
The  initial  distribution  . is  given  by  a sum  of  a bulk  equilibrium  and  a "tail"  Maxwellian;  i.e.. 


f,H-  ^ ^ I T- 


where 


L = "o 


^ 2 

2 

exp 

—m  .V 

~Yf^' 

and 


. 

.1  2 

2 

2 

1 

1277 

1 ■) 

T.tT'u] 

2 7'„r 

>T  = »T 


with  rij-  <3c  n„  and  Tf^j-  » Tj^j-  and  T,,  = 7^^.  We  expect  that  the  dominant  contributions  to  the 
growth  rate  yj^  come  from  2 = 0,  —I  gyroresonances  because  the  i = +1  resonance  is  at  the  nega- 
tive velocity  side  where  the  distribution  function  is  relatively  small.  Solving  the  dispersion  relation  with 
« I and  keeping  only  the  i = 0 and  —I  terms,  we  obtain; 
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where  — (2r,//n,  is  the  perpendicular  ve\oc\iy  of  the  bulk  and  also  the  tail  of  the  distribution 
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function,  and/=  2rr^  <A'xVx/  The  first  and  second  term  of  (3)  correspond  to  Landau  damping  by  the 


1—0  and  1 — —I  terms  respectively.  The  third  term  is  the  driving  term  for  the  instability  and  is 
directly  proportional  to  the  density  of  the  tail.  The  most  unstable  anomalous-Doppler-shifted  mode 
has  kj.  — •'/2k II . 


Let  us  examine  the  circumstances  under  which  these  modes  become  important.  We  can  estimate 
the  resonance  width  in  velocity  space  of  the  I 0 and  Jt  = — I resonances,  given  the  maximum  tail 
velocity  v„.  In  order  for  > 0,  the  third  term  of  (3)  has  to  be,  bigger  than  the  first  term,  i.e.  Lan- 
dau damping  (the  second  term  usually  is  not  as  important).  To  avoid  Landau  damping,  we  usually  need 
> 3v,;,.  But  at  resonance,  v„,  — -I-  <D^)/k^,  so  <«n/k||  < wnv„/((u„  -I-  <u^).  Thus  for 

the  1—0  resonance,  iv,/,  < vn  < a>i^v„/{<o„  •<■««*),  and  for  the  1 — —I  resonance,  we  have 

3v,^  (eo,,  -I-  <«*  )/con  < V||  ^ In  ATC,  for  example  with  =:  and  for  runaway  electron  en- 
ergy with  y„  — lOv,^  and  — ■'/2k, | (which  corresponds  to  maximum  growth),  we  find  for  1 — 0: 

3v,/i  < V||  ^ 3.7v,;,,  and  for  1 — —I  we  find:  8v,^  < v,,  < lOv,^.  For  7",  = IKeV,  this  will 

correspond  to  a runaway  tail  of  about  140  KeV. 


Once  the  anomalous-Doppler-shifted  instability  is  excited,  the  self-consistent  fluctuating  electric 
fields  will  pitch-angle  scatter  the  tail  electrons.  This  process  tends  to  isotropize  the  tail  of  the  runaway 
electron  distribution.  Since  a threshold  exists  for  the  1 — —I  resonance,  there  is  no  pitch-angle 
scattering  outside  the  resonant  region  of  v,, . As  the  high -parallel-energy  particles  are  continuously  be- 
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ing  isotropized,  this  causes  an  increase  of  electron  density  near  the  threshold  of  the  / •>  — | reso- 
nance. As  was  first  pointed  out  by  Papadopoulos this  process  modifies  the  runaway  distribution 
function,  forming  a bump-on-tail  distribution.  This  result  is  confirmed  by  both  particle  simulations  and 
numerical  Fokker-Planck  solutions 


The  growth  rate  y*  and  w * for  the  two-dimensional  bump-on-tail  case  are'*: 
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where  is  the  ion  plasma  frequency,  v*  is  the  velocity  of  the  bump,  and  ReZ'iC,)  is  the  real  part  of 
the  derivative  of  the  Plasma  Dispersion  function.  Depending  on  the  angle  of  propagation,  covers  a 
range  from  to  <o^. 
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Ill  NUMERICAL  MODEL 
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and  y/^  is  given  by  (3)  and  (5). 


(9) 


Eq.  (2)-(9)  are  solved  numerically  in  polar  coordinates  using  an  alternating-direction- 
implicit  method.'*”  Since  it  is  impractical  to  cover  all  directions  of  the  wave  vector  k,  we 
choose  to  solve  the  case  where  is  maximum,  i.e.,  fci  -=  >/2  . We  hasten  to  point  out  that 

this  does  not  mean  we  are  restricting  to  one  mode,  because  for  this  chosen  propagation  angle, 
the  magnitude  of  k still  varies. 

IV  NUMERICAL  RESULTS 


In  previous  work  ” ” we  examined  the  consequences  of  the  bum-on-tail  structure.  Con- 
siderable interest  has  since  been  shown  in  the  actual  formation  of  this  bump-on-tail  distribu- 
tion. We  therefore  present  here  two  types  of  numerical  experiments;  one  which  shows  how 
the  bump-on-tail  distribution  forms,  and  one  which  shows  the  complete  dynamics  of  the  elec- 
tron distribution. 

We  shall  perform  these  calculations  for  Adiabatic  Toroidal  Compressor  (ATC)  experi- 
mental parameters.  Our  theory  predicts  that  the  electrons  of  interest  are  only  mildly  relativistic, 
so  neither  theory  nor  computation  need  be  concerned  with  relativistic  corrections.  The  Alcator 
experiment,  in  contrast,  has  lu,,  » and  the  bump-on-tail  electrons  are  expected  to  be 

highly  relativistic. 

We  begin  by  artificially  removing  the  bump-on-tail  process  from  the  physical  model,  to 
show  how  the  ADS  modes  can  produce  the  bump  on  the  tail  of  the  distribution  function.  We 
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model  the  ATC  experimental  conditions  as  follows:  (o„  “ tUpe.  “ n,  “ 2x10 cm 
riT/fif.  —0.03,  r,  - - 1 kev,  and  - 10.  Here  v„  is  the  maximum  velocity  of  the 

runaway  tail,  and  is  the  runaway  (Dreicer)  electric  field. 

Figure  I shows  the  initial,  monotonically  decreasing  distribution  function,  and  the  one 
into  which  it  evolves  in  10  sec.  Note  the  regions  of  I — 0 and  i — — 1 resonance.  The 
pitch-angle  scattering  due  to  these  resonances  produces  the  bump.  As  we  will  see  below,  the 
bump-on-tail  modes;  when  present,  will  quickly  remove  a bump  of  this  magnitude.  For  now, 
our  purpose  is  to  show  that  such  a bump  can  spontaneously  form  from  a runaway  distribution, 
and  to  examine  its  parametric  dependence. 

In  Fig.  2,  we  change  <Up,, /<*>,.,,  = 1 to  = 1.25  which  corresponds  to  higher  densi- 

ty. We  see  that  the  threshold  (the  bump)  gets  closer  to  the  thermal  region.  Eventually  as  the 
density  increases  the  bump  may  get  so  close  that  it  may  merge  with  the  thermal  region,  and 
the  bump-on-tai)  modes  could  no  longer  be  excited. 

It  is  the  flattening  of  the  bump,  and  the  resulting  bump-on-tail  modes,  which  accounts 
for  the  slowing  down  of  the  runaway  electrons.  Therefore,  we  expect  that,  in  high  density 
operations  where  the  bump  approaches  the  thermal  distribution  the  runaway  electron  energy  is 
higher,  since  the  bump-on-lail  modes  provide  less  energy  loss. 

When  we  switch  on  the  bump-on-tail  mode,  we  do  not  expect  to  see  the  bump  anymore, 
because  the  bump-on-tail  mode  occurs  on  a much  faster  time  scale  and  flattens  any  bump.  In 
Fig.  2,  we  show  the  time  development  of  ,/  (v,,  ) with  both  the  anomalous-Doppler-shifted  and 
bump-on-tail  modes  present.  The  parameters  are  the  same  as  Fig.  1 except  nf/n„  = 2%.  The 
dotted  curve  shows  the  distribution  function  at  steady  state,  about  I /asec.  after  initiation. 

Note  that  no  bump  is  formed,  and  /(v,,  ) is  flatter  at  the  tail  than  initially.  We  find  the 
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anomalous-Doppler-shifted  instability  dies  down  because  of  the  depletion  of  particles  by  the 
1 - -1  resonance.  This  is  consistent  with  (3)  which  says  the  growth  rate  is  directly  propor- 
tional to  the  magnitude  of  /at  the  i “ — I resonance. 

We  have  also  calculated  the  synchrotron  radiation  power  due  to  uncorrelated  electrons 
and  find  that  the  power  increases  by  about  75%  in  about  1 jisec  (See  Fig.  4)  and  levels  off  after 
the  instabilities  die  down.  This  is  consistent  with  experimental  observations.^  We  conclude 
that  the  excitation  of  the  anomalous-Doppler-shifted  mode  is  definitely  related  to  the  increase 
of  synchrotron  radiation.  If  we  follow  /(v,,  ) longer,  we  expect  the  resonance  region  of 
£ = -1  to  be  filled  up  again,  because  E,,  can  accelerate  particles  from  low  v,,  to  high  v,,. 
When  this  happens,  the  instability  should  be  excited  again  and  the  synchrotron  radiation 
should  also  increase.  We  have  followed /(v„  ) up  to  2.8  x 10  sec  when  the  tail  region  is  al- 
most filled  up  again  (see  the  circles  of  Fig.  2).  We  find  the  instability  is  not  excited  because  the 
negative  slope  at  f = o also  increases  tremendously.  At  first  we  believed  this  was  due  to  fi, 
pushing  the  bulk  electrons,  but  we  found  this  effect  still  remained  after  we  switched  off  ^n. 
We  conclude  this  must  be  due  to  isotropization.  In  fact  it  can  be  easily  estimated  that  complete 
isotropization  of /(v|,  ) at  v,)  = 3 takes  about  1 msec.  Therefore  in  about  a third  of  this 
time,  the  slope  can  increase  enough  to  damp  out  the  anomalous-Doppler-shifted  mode.  Exper- 
imentally it  is  observed  that  the  synchrotron  radiation  increases  in  steps  a few  msec  apart.  If 
the  increase  in  radiation  is  related  to  the  excitation  of  anomalous-Doppler-shifted  mode  which 
we  think  it  is,  then  v„  and  the  instability  threshold  for  anomalous-Doppler-shifted  mode  must 
increase  so  that  the  -=  0 resonant  region  can  reach  4 v,^  or  5 v,/,. 

V SUMMARY 

We  have  developed  the  analytic  tools  for  describing  the  effect  of  anomalous  Doppler- 
shifted  plasma  waves  and  bump-on-tail  modes  upon  the  tokamak  electron  distribution  func- 
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tion.  These  effects  have  been  self-consislently  incorporated  in  a numerical  solution  of  the 
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Fokker-Planck  equation  as  a quasi-linear  diffusion  ope/ator.  Based  on  results  using  this  model, 
we  find  the  following  general  properties  of  the  tokamak  electron  distribution: 

1.  The  onset  of  the  anomalous-Doppler-shifted  modes  gives  rise  both  to  a suprath- 
ermal  synchrotron  radiation  spectrum  and  formation  of  a bump-on-tail  distri- 
bution. 

2.  The  combination  of  anomalous-Doppler-shifted  and  bump-on-tail  modes  are 
necessary  to  slow  down  the  runaway  electrons. 

When  comparing  different  experiments,  our  studies  predict: 

3.  Experiments  at  higher  electron  number  density  should  show  lower-amplitude 
anomalous-Doppler-shifted  modes,  and  higher  mean  runaway  energy. 

4.  At  any  density,  the  threshold  in  v for  onset  of  the  anomalous-Doppler-shifted 
mode  shifts  outward  with  time. 

Finally,  the  sudden  reduction  in  Vn  caused  by  onset  of  the  bump-on-tail  modes  and 
resultant  flattening  of  the  electron  distribution  function  may  cause  voltage  spikes.  This  loss  in 
V||  takes  place  with  minimal  loss  of  energy,  breeding  very  hot  electrons  with  large  vx  and 
small  V|| . Such  electrons  would  be  trapped  in  poloidal  magnetic  field  minima,  and  could  drift 
to  the  outside  of  the  system,  leading  to  liner  damage  such  as  seen  in  TFR  and  Alcator. 
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V||/V*, 

Figure  3 — /(V||  ) at  / — l after  both  anomalous- Doppler-shifted  mode  and  bump- 
on-tail  mode  die  (dotted  curve),  /(vn  ) at  r — 2.8  x 10  sec  (circle  curve), /(v,,  at  r —0  is 
given  by  solid  line,  (i»„/a»p,  “I,  * 2 x lo*^  cm  T,  — I keV,  nj/rtg  — 2%,  — 

0.06  and 
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